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Abstract 

In the random case setting, scrambled polynomial lattice rules as 
discussed in p] enjoy more favourable strong tractablility properties 
than scrambled digital nets. This short note discusses the application 
of scrambled polynomial lattice rules to infinite-dimensional integra- 
tion. In [5], infinite-dimensional integration in the random case setting 
was examined in detail, and results based on scrambled digital nets 
were presented. Exploiting these improved strong tractability prop- 
erties of scrambled polynomial lattice rules and making use of the 
analysis presented in [5] , we improve on the results that were achieved 
using scrambled digital nets. 

1 Introduction 

In a recent series of papers, [2], [9], [5], [11], the problem of infinite-dimensional 
quadrature has been studied. Such problems have many applications, e.g. in 
mathematical finance, see [10], [8], where series expansions are used to rep- 
resent particular random variables. 

This short note focuses on the random case setting, which was addressed 
in [5]. In the latter paper, a complete and general analysis was presented, 
clearly showing the reader how to employ a given quadrature rule. Further- 
more, results based on the scrambled Niederreiter sequence were presented. 
Recently, [T], it was shown that for multivariate integration in the random 
case setting, scrambled polynomial lattice rules possess more favorable strong 
tractability properties than any scrambled digital sequences. 
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This raises the natural question, whether the results based on the scram- 
bled Niederreiter sequence can be improved using scrambled polynomial lat- 
tice rules. In this short note, we give an affirmative answer to this question. 
In particular, the contribution of this note is the following: Using the anal- 
ysis presented in [5], we employ the multivariate in tegration result from jl] 
to improve on the results presented in [3]. Furthermore, due to the growth 
of the t-value of order t x s, it is clear that the results as presented in this 
paper cannot be achieved using digital sequences. Finally, we remark that 
results on multivariate integration using lattice rules in the random setting 
have not appeared in the literature. 

We place ourselves in the same setting as discussed in [5J, use the same 
algorithms, but employ a different quadrature rule, in particular with better 
strong tractability properties. In the interest of giving due credit to the 
authors of [5] and also of saving space, we have decided to proceed as follows: 
We very briefly recall the function space introduced in [5] and the sampling 
regimes, but introduce cost and worst-case errors under the premise that 
algorithms are based on scrambled polynomial lattice rules. The interested 
reader is referred to [S] for a complete and general treatment of the problem 
studied in this note. 

2 The Setting 

In this section, we briefly recall the function space and the sampling regimes, 
cost and errors as introduced in [5]. Regarding notation, as in [5], v is used 
to denote finite subsets of N, the set {1, . . . , s} is denoted by 1 : s, lastly, we 
write Xk di Uk for sequences of positive real numbers Xk and yk, if Xk < cyk is 
valid for k G N and some constant c > 0. 

2.1 The Function Space 

We briefly remind the reader how to construct functions of infinitely many 
variables, as presented in [5]. Essentially, we start with a one-dimensional 
reproducing kernel Hilbert space, construct spaces of finitely many variables 
as tensor product spaces and take limits to allow for infinitely many variables. 
Coordinate weights j v = Y]_j ev Jj, v ^> which indicate teh importance of 
the variables Xj, j 6 v, ensure convergence of the relevant quantities. 



In particular, we consider the reproducing kernel 

H x > V) = 3 + ( x2 + V 2 )/ 2 ~ max(x, y) , 

x, y G [0, 1], and consider the Hilbert space H(l + jk), for a weight 7 > 0, 
whose norm satisfies 

1 f(y)dy] +7- 1 f\f') 2 (y)dy. 

/ Jo 

To allow for functions of finitely many variables, we consider the reproducing 
kernel 

K v {x, y) = + 7jk(xj, yj)) 
and of course the associated Hilbert space H(K V ) is of tensor product form 



H(K v ) = (g)H(l + lj k). (1) 

To define functions of infinitely many variables,w e define the measurable 
kernel K on [0, 1] N x [0, if 

K(x,y) = ^2>y v K v (x,y) = k(xj, yj) , 

for x, y G [0, 1] N and denote the associated space by H(K), which, see [5j 
Lemma 6], consists of all functions 

/ ^ ^ fvi fv G H v , 

V 

for which 

T,^ 1 \\M\l<oc, 

v 

and, in case of convergence, 
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2.2 Sampling regimes, cost, and worst-case error 

In this subsection, we introduce randomized algorithms for the integration of 
functions / : [0, 1] N — > R; the reader is referred to [15] and [2] for a detailed 
discussion. 

Follwoing [2], [5], two sampling regimes, which specify the domains from 
which the integration nodes can be specified, are introduced. 

Fixed subspace sampling restricts this domain to a finite-dimensional affine 
subspace 

X»,a = {x E [0, 1] N : Xj = a for j G N \ v} 

for a finite set I^i/CN and a G [0, 1]. We remind teh reader that essentially 
one only specifies those co-ordinates included in v, the remaining ones are 
specified via teh anchor point a. 

Variable subspace sampling generalizes this idea to a sequence of finite- 
dimensional affine subspaces 

X Vl: a (— ^-1)2)0, C . . . , 

where v = (ui)teN is a given increasing sequence, Vi C N, and a G [0, 1]. This 
sampling schemes allows us to choose integration nodes from subspaces of 
different dimensionality. To be able to compare introduce a cost function c, 
which si a mapping 

c : [0,1] N ^ NH {oo} 

and which is to quantify how costly it is to evaluate teh integral / at teh 
integration nodes. Following [2], [5], we formulate the cost of evaluating / 
at the integration node x in terms of the dimension of teh finite-dimensional 
subspace from which the integration node is chosen, thsi means for fixed 
subspace sampling 

f \v\ , if x G X v , a , 2 \ 
v ' a \ oo , otherwise. 

For variable subspace sampling, which allows for the integration nodes to 
be chosen from a sequence of finite-dimensional subspaces, we choose the 
subspace with the smallest dimension in which the node lies, 

c v ,a( x ) = inf {dim(A;. )a ) : x G X ViA } , (3) 

and set inf = oo. 
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The randomized quadrature formulas employed in thsi note are based on 
scrambled polynomial lattice rules, 

b m 
i=l 

where X{ £ [0, l) s are obtained by scrambling a polynomial lattice rule, see 
PP. Defining 

(*v,af)(x) = f(xv,a), (5) 
we denote teh randomized quadrature formulas of interest in thsi note by 

Qn,s,a = Q[log b (n)],b,l:s ° ^l:s,a , (6) 

where n denotes the number of points of which the quadrature rule is com- 
prised. Intuitively speaking, we carefully specify dimensions 1 to s via scram- 
bled polynomial lattice rules and employ the anchor a £ [0, 1] for the sub- 
sequent dimensions. It is clear from the previous discussion, that for fixed 
subspace sampling, the notation introduced in ([6]) is sufficient. For variable 
subspace sampling, however, we allow our integration nodes to be chosen 
from subspaces of different dimensions, in this sense the letter s is not suffi- 
cient. We remark that this is addressed in Section |H 

Next, we wish to discuss the cost of the randomized algorithms. As we 
only discuss randomized algorithms based on scrambled polynomial lattice 
rules in this note, we define the cost of fixed and variable subspace sampling 
under the premise that the randomized algorithm is based on a scrambled 
polynomial lattice rule. This simplifies the discussion, the cost model em- 
ployed in [5], which stems from [2], allows for a much more general class of 
algorithms, see also [9] for an even more general cost model. 

Essentially, the cost of evaluating a randomized algorithm Q is given by 
the sum of the costs of evaluating the function at the integration nodes cho- 
sen from the finite-dimensional subspaces. For teh fixed subspace sampling, 
assuming that nodes are chosen from Xi :s>a 

cost f ix (Q) = ns . 

For variable subspace sampling, we choose our integration nodes from 
a sequence of finite-dimensional affine subspaces, index v = (fj)™^, where 
m < n, as we use an n-point quadratuer rule, where Vi — 1 : s^, i = 1, . . . , m. 
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For the subspace indexed vy eth integration nodes would be based on 
scrambled polynomial lattic erules whose integration nodes lie in [0, l) Si , and 
we denote the number of those integration nodes by n v ., where of course 
JZiLi n vi = n. Consequently, we have 

m 

COSt var (Q) = ^ S i U v l ■ ( 7 ) 
i=l 

Finally, following [5], for class of integrands F, where / G F is a mapping 
/ : X R, X C [0, 1] N , we use the notation 



1(f) = [ f(x)dx 

J X 



and denote the worst-case error of a randomized algorithm Q, used to ap- 
proximate integrands / in the class F by 

e(Q,F) = sup (E (1(f)- Q(f)f) 1/2 . 
feF 

Lastly, minimal errors, which are of great importance in information- 
based complexity, [15], [12], [H], are defined by 



e NJix (F) = inf {e(Q,F) : cost fix (Q,F) < N} 

and 

(F) = inf {e(Q, F) : cost var (Q, F) < N} . 

The following result on numericla integration in H(Ki :s ), see Equation 
([I]), stems from [1]. 

i 

Theorem 2.1 Assume YlJLilj'" < °°> fore > 0, then 
e(Q b , m ,i:s,H(K 1:a )) < c £ n-^ 2 ~ e \ 

where n = b m and Qb,m,i-.s is a scrambled polynomial lattice rule as defined 
in Equation (EH). 

Proof. From the proof of [T6J lemma 7], it is clear that the function space 
H(K Vs ) can be embedded in the space Vi iSj7 , as defined in [TJ, from which 
the result follows immediately. □ 

Concluding this section, we point out, that the forthcoming results will 
be presented for functions in B(K)), the unit ball in H(K). 
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3 Results on Fixed Subspace Sampling 

To fully specify the fixed subspace sampling algorithm, we only need to spec- 
ify the dimension of the finite-dimensional subspace employed for sampling, 
and the number of integration nodes, which are based on a scrambled poly- 
nomial lattice rule, that we employ. As we wish to minimize worst-case errors 
for a fixed bound on the cost, say N, both, the dimension and the number 
of integration nodes are functions of N. 

Corollary 3.1 (Corollary 1, [5]) Let e > and let jj x j~ a , a > 3. 
Choose 

a — l 

n x jV"+2-E 

and 

3-E 

s x N a + 2 -c 
for N 6 N. Then, for Q N = Q n ^ a 

(3-£)/2(a-l) 

e{Q N ,B{K))<N 

and 

cost fix (Q N ,B(K)) . 

Proof. The result follows immediately from [U Theorem 1], where we set 
a > 3. □ 

Remark 3.1 In the same result was established under teh stronger as- 
sumption on the weights jj x j~ a , a > 4. The result presented in Corollary 
\3.1\ is optimal for a > 3, see JE, Corollary 3}. 

4 Results on Variable Subspace Sampling 

We carry out variable subspace sampling using the so-called multi-level ap- 
proach, which was first introduced in [3], pE], see also [6], [7]. teh idea under- 
lying the multi-level approach is the following: We fix a sequence of sets 

v x C ■ • • C v L 

and the associated finite-dimensional affien subspaces 
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We use the integral associated with the finite-dimensional subspace of the 
largest dimension, I(^ VLA f) to approximate /(/). However, we rewrite 
I(^ VL>a f) as follows 

L 
1=1 

setting \l/„ 0i a/ = 0- Each of the integrals I(^ Vl>a f — W Vl _ lia f) is now approx- 
imated using an independent randomized algorithm, based on a scrambled 
polynomial lattice rule, in particular we use a randomized algorithm 

L 

o(/) = X)Q»«.'«.»(/-*^-i.*/)' ( 8 ) 
i=i 

so at level I, we use an algorithm based on a scrambled polynomial lattice 
rule consisting of &L lo g&( n i)J points, which lie in [0, l) s '. The error associated 
with this algorithm can be split into bias and variance, 

E (1(f) - QU)f = (1(f) - I(*i:s L ,af)) 2 + Var(Q(/)) , 
in particular 

L 

Var(Q(/)) = Var {Q ni , si>a (f - *i: Sl _ ua f)) , 
i=i 

see [5]. Regarding the cost, from Equation (j7|), 

L 

COSt var {Q, B(K)) = ^ S l n l ■ 
1=1 

By definition of variable subspace sampling, the dimension si increases with /, 
but one would expect the variances Var (Q ni)Sl) a(f — ^i:a } _i,o/)) to decrease 
as / increases; the challenge is to trade of these effects well. 

Corollary 4.1 (Corollary 4, |5j) Assume that jj x j~ a , for a > 3, let 
< e < min(6, a — 3) and put 
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Choose L, si, ni according to Equations (26), (27), and (28) in |3J/ ; and 
let a € [0, 1]. Take the corresponding multi-level algorithm Qn according to 
Equation (JS} based on the scrambled polynomial lattice rule. Then 



e(Q N ,B(K)) ± 

and 



7V-( 3 - £ )/ 2 , if a > 10, 

^-(3^/23=1 ) if a < IQ , 



cost var (Q N ,B(K)) ■< N . 

Proof. The proof follows immediately from [3 Theorem 4], with a 1 = 3 + e. 

□ 

Remark 4.1 The same error bounds were established in JSjj, but the rate 
was only established for a > 11, whereas here it is achieved for 
a > 10, due to an improved strong tractability result. Of course, for a > 10, 
this result is optimal. Furthermore, we conclude that for (at least) a > 7, 
variable subspace sampling improves on fixed subspace sampling. 

Remark 4.2 We alert the reader to ITI^ . where infinite-dimensional inte- 
gration in the worst-case setting is studied. In [77]/ . rank-1 lattice rules are 
employed as a basis for the algorithms, and we remark that in the worst-case 
setting, polynomial lattice rules have not been shown to improve on rank-1 
lattice rules. 



5 Future Research and Applications 

Applications of infinite-dimensional quadrature problems abound, to name 
but two from mathematical finance that received attention at MCQMC 2010, 
we list series expansions, such as the Karhunen-Loeve expansion in the case 
of Brownian motion, see e.g. [10] for a related publication, and the shot-noise 
representation in the context of Levy processes, see [Hj . 

In the case of infinite-dimensional quadrature, the theory, as presented 
in [2], [9], [5], [11], is ahead of the applications in the following sense: If 
we can specify the decay of the weights, all parameters of our algorithm are 
readily determined. However, in practice, the study of the weights has still 
not received enough attention, but it is clear that in order to fully reap the 
benefits of the algorithm, this needs to be addressed. 
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